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O ■ Abstract 

In [I], Baumann et al. presented a new formalism for studying cosmological 
systems where the characteristic scale of non-linearities, fcw L , is much smaller 
than the Hubble scale l~L~ l . By integrating out the short-wavelength modes, 

£T) ' it is possible to obtain an effective theory of long-wavelength perturbations 

that is described by an imperfect fluid evolving in an FRW background. As 

t— ( ■ the long-wavelength perturbations remain small even when the short-scale dy- 

0^ ' namics are non-linear, the tools of linear perturbation theory may be applied. 

The work in pQ deals only with matter in the form of a pressureless perfect 
fluid with zero anisotropic stress, and also assumes that the short-scale gravi- 

. !_< | tational dynamics are Newtonian. In this work we extend this formalism to the 

case of a perfect fluid with pressure, and in particular to the case of preheat- 
ing after inflation, where the matter content of the universe can be modeled 
by two coupled scalar fields. We discard the assumption that the short-scale 
gravitational dynamics are Newtonian. We find that our results differ from 
Baumann et al.'s even when the pressure is set to zero, which suggests that 
relaxing their assumptions creates appreciable changes in the long-wavelength 
effective theory. We derive equations of motion for the total density pertur- 
bation and matter velocities during preheating, as well as linearized Einstein 
equations for the long-wavelength metric perturbations. We also present the 
equations governing the effective long- wavelength scalar field dynamics. 
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1 Introduction 



Recently, Baumann et al. have proposed a formalism for analytically studying long- 
wavelength dynamics even when the density contrast 5 of a universe grows large pQ. 
Their method applies when the characteristic scale of non-linearities, k^, is much 
smaller than the Hubble scale "H -1 . This hierarchy allows us to integrate out the 
short-wavelength modes by smoothing all perturbations over a scale A -1 that lies 
between the non-linear scale and the scale of the long- wavelength perturbations. This 
procedure gives an effective theory of long-wavelength perturbations that is modeled 
to lowest order by an imperfect fluid with effective energy-momentum pseudotensor 
Tfj_ u , evolving in an FRW universe. The properties of the effective fluid are determined 
by the interactions of the short- wavelength modes. 

The equations of motion in the effective theory can be expressed as a derivative 
expansion of the long-wavelength variables, with higher order terms being suppressed 
by powers of (k/A) 2 or (k/k NL ) 2 . Since the long-wavelength perturbations during 
preheating remain small even when the underlying dynamics have become non-linear, 
we can then apply perturbation theory techniques. The coefficients of the effective 
theory can either be obtained by matching to the results of numerical simulations, or 
left as free parameters to be matched to measurements from experiment. 

Baumann et al. considered a universe filled with cold dark matter, modeled as 
a pressureless perfect fluid. In this paper we consider a universe where the matter 
content is a perfect fluid with pressure: in particular, we consider matter in the 
form of two coupled scalar fields. We also make fewer assumptions than Baumann 
et al. For example, we discard their assumption that the short-scale gravitational 
dynamics are Newtonian. We obtain somewhat different results from those given in 
PQ, even when we set the pressure to zero in our equations, suggesting that relaxing 
the underlying assumptions on the short-scale gravitational dynamics significantly 
affects the long-wavelength perturbations. We hope that our approach will thus be 
applicable in more general scenarios than a ACDM universe. 

One particular example where our work applies is the period of preheating fol- 
lowing inflation. Inflation ends when the slow-roll conditions are violated, and the 
inflaton begins to oscillate around its ground state. As the inflaton field is coupled 
to the standard model (SM) matter fields, these oscillations cause the energy stored 
in cf) to be transferred to the SM fields. This process is known as reheating^. 

Reheating was originally analyzed using perturbative quantum field theory, until 
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it was realized that the coherent nature of the inflaton field at the end of inflation 
renders this picture inaccurate. The excitation of SM fields during reheating was then 
reformulated as a semi-classical problem, in which the quantum mechanical produc- 
tion of SM matter particles takes place in the classical background of the inflaton 
field. The most common toy model used to describe this situation couples a scalar 
matter field \ to the inflaton field <fi. In this model of reheating the \ field is excited 
via parametric resonance, and the inflaton decays rapidly as a result, so that the 
dynamics quickly become non-linear. In this regime the standard tools of linearized 
cosmological perturbation theory no longer apply, so most studies of reheating have 
relied on numerical simulations [£]. In particular, the occupation numbers of x an d 4> 
are found to be large in the non-linear domain, so that they may be treated as classi- 
cal fields and studied using lattice simulations. This period of rapid energy transfer 
from <fr to x is known as preheating. It results in a highly non-thermal distribution of 
matter fields that then thermalizes to give the initial conditions for a hot big bang. 

Although computer simulations are useful for providing visual representations of 
preheating dynamics, it would be desirable to have an analytical description of the 
non-linear dynamics. In addition to allowing us to make explicit calculations, it could 
also offer a greater understanding of the scalar field interactions. As preheating is a 
highly inhomogeneous process, the density contrast 5 = p/p — 1 for both <fi and x soon 
exceeds 1, so that conventional cosmological perturbation theory (in which we expand 
in 5, matter velocity v, and metric perturbations $) cannot be applied. However, 
numerical studies have shown that there is a large hierarchy between the non-linear 
scale and the Hubble scale during preheating, with ~ 0.01H _1 [4]. Moreover, the 
metric and velocity perturbations remain small, even while the inhomogeneous parts 
of and x grow large [5]. Thus it seems that the period of preheating provides a 
suitable test case for setting up an effective theory of long-wavelength perturbations. 

This paper is organized as follows. In Section [2] we describe the rules for expanding 
the equations of motion in our formalism. In Section [3] we give the basic equations 
that govern the scalar field and metric dynamics during preheating. In Section H] we 
derive the effective theory of the long- wavelength modes, and in Section [5] we give the 
evolution equations for the long- wavelength theory. We conclude in Section [6j 
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2 Cosmological Perturbation Theory: A Velocity 
Expansion 



In cosmological perturbation theory, we usually expand to linear order in the density 
perturbation Sp, matter velocity v, and metric perturbations $, so that linear per- 
turbation theory is no longer valid when the density contrast 5 > 1. In our work, we 
follow the alternative approach outlined in pQ, where we expand up to order v 2 in 
velocity and metric perturbations, and we do not expand in Sp. 

In perturbation theory, we find that v is related to the Newtonian potential $ 
and the density contrast 5 by 

v 2 ~ <M (2.1) 

In the non-linear regime, when 5 ~ 1, we therefore find that v 2 ~ $. 

On small scales, gradients of the gravitational potential can change the power- 
counting of standard perturbation theory, as the short-wavelength modes have large 
momenta k. The net result is that, at the non-linear scale, each gradient of $ reduces 
the order in v by one. Therefore, when we carry out an expansion to order v 2 , we 
expand to linear order in $ and second order in V$. 

Simulations of the scalar fields during preheating have shown that the metric 
fluctuations and the gradient energies of the fields remain small, even when the fields 
themselves (and their kinetic energies) grow large [5]. Therefore, it is valid to apply 
this expansion when studying preheating dynamics. 



3 Preheating: The Basic Equations 

We will consider one of the simplest models of preheating, in which a massive inflaton 
field 4> with inflaton potential V\ = |m 2 2 couples to a massless scalar field x through 
a potential V2 = \g 2 <t> 2 X 2 - (Our approach can easily be extended to other models with 
different interaction potentials.) The full potential is therefore 

V(<j>, X ) = V 1 + V 2 = Kn 2 4> 2 + ^ W (3.1) 
At the end of inflation is a homogeneous field that approaches 

(/) -» p sin rot, (3.2) 
V37rmt 
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where t is the cosmic time. The spacetime is described by a perturbed FRW metric 

ds 2 = a 2 ( V ) (-e 2 W + e" 2 *dx 2 ) , (3.3) 

where rj is the conformal time, we have ignored vector and tensor perturbations, and 
we use the Poisson gauge. In the case of zero anisotropic stress, we find that \& = $ 
to first order. 

The basic equations that describe the dynamics during preheating are the Klein- 
Gordon (KG) equations for and and the Einstein equations. The KG equations 
are: 

□0-^ = (3.4) 
n x -V x = (3.5) 

where = ^r. The energy- momentum tensor of the scalar fields is given by 

T» = d^d^ + d^X ~ K Q^0«9> + l -d aX d a X + V) (3.6) 
Furthermore, T£ obeys the conservation law 

= 0. (3.7) 
Finally, the Einstein equations are given by 

= 8ttGT». (3.8) 

In the rest of this section we expand the KG equations and the continuity equations 
to order v 2 in the metric perturbations and the matter velocities, but we do not expand 
in p or P. We find that the Einstein equations have a very similar form to those given 
in PP, at least to leading order in long- wavelength perturbations. We therefore give 
these equations in Appendix [Bj In order to apply the expansion to the system of 
coupled scalar fields described above, we first consider the ways in which this system 
may be interpreted as a sum of perfect fluids. 

3.1 The Scalar Field as a Perfect Fluid 

It is possible to interpret a system of interacting scalar fields as a system of interacting 
fluids [7J EJ HI [10] . We can therefore expand the equations of motion in terms of p 
and v, or in terms of and x- This will give us two equivalent descriptions of the 
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matter fields during preheating, and we can go back and forth between the two using 
the correspondence outlined below. It will be useful to consider both descriptions, as 
they offer different physical insights and different calculational advantages. The fluid 
description allows us to consider the energy density of the entire system using only 
one variable, p, and gives a simpler formulation of the equations of motion. The scalar 
field description allows us to consider the individual field perturbations directly. 

The energy momentum tensor of a perfect fluid with energy density p, pressure 
P, and zero anisotropic stress is given by 

T? = (p + P)u»u u + 6£P, (3.9) 

where w M is the instantaneous 4- velocity of the fluid. The components of the 4- velocity 
are related to the matter velocity v by the equations 

u° = a- 1 e-* 7 (w), u< = crW (3.10) 

Mo = — ae^'y(v), m = ae _4 V, (3-H) 

where 7(f) := (1 — v 2 )^ 1 ^ 2 . Comparing Eq. (13.91) with the form of in Eq. (l3.6|) . we 
see that the interacting scalar fields d> and x can be treated as the sum of two "kinetic 
fluids" with energy density and pressure [8] 

p^ := -~cW>c> a Px :=-l daX& x x (3.12) 

P^:=~d a ^ P x -=-\d aX d a X (3.13) 

and a single "potential fluid" with energy density and pressure 

pv := V (3.14) 
P v := -V. (3.15) 



-d»d> 



<:=^fc= «5 = ^4= (3.16) 



and to order v 2 the matter velocity v is given by 

&6 did 



d°cj) d <p x 
The 4-velocity is not defined for the potential fluid. 



iven 


by 




-d^x 


v- 


■d a xd a x 


d l x 


diX 


d°X 


doX 



v* = = (3.17) 
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The total energy-momentum tensor is given by the sum of the individual energy- 
momentum tensors of these three fluids. The total veocity perturbation is given by 

v = — (3.18) 

where p := p^ + p x + pv and P := P^ + P x + Py are the total energy density and 
total pressure respectively. 

The physical interpretation of the scalar field as a fluid is valid as long as 
is timelike. In our case, we are assuming that v is small compared to 5, which is 
equivalent to assuming that the gradient energy of the scalar fields is small compared 
to their kinetic energy. Thus it is reasonable to assume that remains timelike during 
preheating. The form of shows that the anisotropic stress is zero. Therefore, we 
can set the metric perturbations $ = \1/ to first order, and will do so for the rest of 
the paper. 

3.2 The Energy-Momentum Conservation Equations 

The components of the energy-momentum tensor calculated to order v 2 in terms of 
p, v are: 

T ° = -(p + P) 7 2 + P (3.19) 
Tq = -v\p + P) (3.20) 
Tj = (p + P)v l v j + 5)P (3.21) 

In terms of 0, Xi the components are: 

- V (3.22) 
1$ = ^(d^t + dixdox) (3-23) 

T; = ±(d i <f>d j( !> + d i xd j x)-S}v 



In [TJ, the equations of motion of the long- wavelength perturbations are derived us- 
ing the Euler equations in the Newtonian approximation. We will take a different 
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approach and use the full general- relativistic conservation equation (I3.7p . We can 
project this equation along u u , or along the orthogonal direction: 

u v V^ = -~?=^ (V=9(P + p )0 + ^P = (3-25) 

[f v + u u u u )V^ = (p + P)«"V/ + <9 ff P + u u u»d^P = (3.26) 
To order t> 2 , Eq. (13.251) is: 

p + $((/0 + P)*/) - v%P + {?>% - 3$ + 2v • v)(p + P) (3.27) 
To order t> 2 , the spatial components of Eq. (l3.26p are: 

(p + P) (ui + v j djVi + di%) + diP + Vid P + v iV j djP = (3.28) 

3.3 The Einstein Equations 

To write down the Einstein equations to order v 2 , we decompose the Einstein tensor 
G^v into a homogeneous background G^, a part linear in the perturbations 
and part non-linear in the perturbations (G^ L ). We can then write the non-linear 
Einstein equations in the form 

G%, = 8irG(r^ - TfJ) (3.29) 
where the energy-momentum pseudotensor t^ v is given by 

qNL 

V = T, V --^ (3.30) 
The velocity expansion of the Einstein equations is given in Appendix IB. 11 

3.4 The Klein-Gordon Equations 

Expanding the Klein-Gordon equations for 0to order v 2 gives: 

(1 - 2^)<9 O 2 - <9 2 + (2H - - 3$ - ^)<9 O + <9,($ - tydrf + = 0, (3.31) 
with an analogous equation holding for \- 
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4 Integrating Out Short- Wavelength Modes 



In this section we describe how to integrate out the short-wavelength modes to ob- 
tain a long-wavelength effective theory. Integrating out the short-wavelength modes 
amounts to averaging perturbations over a smoothing scale A -1 . Since we are inter- 
ested in the theory at scales k' 1 much larger than the non-linear scale k^, we choose 
a smoothing scale A -1 >> k^. 

The smoothing of perturbations corresponds to a convolution of all fields X 
{ />. $, pv} with a window function W A . We define the long-wavelength mode Xi of a 
field X to be 

X l ee [X] A (x) = j dV W A (|x - xl)X(x'). (4.1) 
The short-wavelength mode X s of X is then defined by 

X = X ; + X S . (4.2) 

We will assume W A to be Gaussian for convenience. We also assume that W A satisfies 
the following conditions: 



d f W A = —djW\ = A 2 (x - x') j W A (4.3) 
v d f W A = didjWA = -A 2 8ijW A 



d v d r W A = didjW A = -K 2 8ijW A + A 4 (x - x')' ; (x - x.') j W A (4.4) 



If we smooth general bilinear and trilinear quantities using W A , we obtain the follow- 
ing results: 

[f9]A = fi9i + [fs9s]A + ^V/, • V 9l + . . . (4.5) 
[fgh\\ = fwihi + [fs9sh s ] A + fi[g s h s ] A + gi[f s h s ] A + hi[f s g s ) A 
+ -j^UNgi ■ + giVfi ■ Vh t + h t Vfi ■ Vgi) 
+ -j^iyfi ■ V[g s h s ] A + Vgi ■ V[/A]a + V/ij • V[f s g s ) A ) 
+ ^2(V 2 fi\g,h,] A + V 2 gi[f s h s ] A + V 2 h[f s g s ] A ) + ... (4.6) 

The explicit calculations used to derive Eq. (I4.5ti4.6p are given in Appendix |A] The 
expressions are given up to higher derivative terms of order k 2 /A 2 , where k is a 
characteristic frequency of a long-wavelength mode. The higher derivative terms are 
therefore suppressed, as k « A. As mentioned in the Introduction, our expression 
for smoothed trilinear quantities differs from that given in pQ, due to the presence 
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of extra terms that are absent in their formula. It is possible that the discrepancy 
is due to extra (unspecified) assumptions that have been imposed in their paper to 
allow these terms to be dropped. 



5 The Long- Wavelength Effective Theory 

We now use the equations from Section [3] and Section 0] to derive the long- wavelength 
effective theory of preheating dynamics. Smoothing the energy-momentum pseu- 
dotensor in f !3.30j) . we find that the effective energy-momentum pseudotensor [t^ u ]\ 
has the form 



where r' depends only on the long-wavelength perturbations, depends on the 
short-wavelength modes, and r^ v contains higher- derivative corrections that are sup- 
pressed by powers of k 2 /A 2 . Throughout this work, we will drop all such higher- 
derivative corrections that result from smoothing. We also drop non-linear metric 
contributions to tL,. 

The pseudotensor [r^ u }\ describes an imperfect fluid. Thus the effective theory we 
obtain after smoothing is an imperfect fluid with density perturbation 5i and velocity 
perturbation v/, evolving in a background FRW metric with scalar metric perturba- 
tions $/ and i. As the scale of non-linearities is much smaller than the smoothing 
scale, these long-wavelength perturbations remain small even when the small-scale dy- 
namics have become non-linear. Therefore, we may apply linear perturbation theory 
in the effective theory even when 5 » 1 at small scales. 

To leading order in the long-wavelength perturbations, we may write the energy- 
momentum pseudotensor as: 



where 6 = d{Ui, ( is the bulk viscosity, and is the anisotropic stress. The pseu- 
dotensor is given as a derivative expansion, with higher order terms being suppressed 
by (h/kpti,) 2 . All of the quantities in (15.21) are quantities in the effective theory: for 
example, p := p c s = p e g + 5p c ff- In order to avoid cluttering the notation, we will 
omit the "eff " subscript from relevant quantities, with the exception of Section 15.11 
when we discuss the renormalization of background pressure, energy density, and 
anisotropic stress due to short-wavelength dynamics. Although we began with a sys- 
tem of coupled scalar fields that had zero anisotropic stress and zero viscosity, we 






Ma = (p + P)u„u v + {p- ce)g^u + £ 
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will find that both anisotropic stress and viscosity are induced in the long-wavelength 
effective theory. 

We use the following ansatz for the anisotropic stress: 

= —r)(Tij, <7ij := V(ij) — -SijVk,k, (5.3) 

where rj is the shear viscosity. To lowest order, the pressure perturbation 5P is 

5P = c 2 sP 5, (5.4) 

where c 2 s is the speed of sound squared. We also define the equation of state parameter 
w = — , and the dimensionless paramter c 2 is that characterizes the viscosity by 

The anisotropic stress induced in the effective theory satisfies 

" ~p~W j = ~ Cms u (5,6) 

To leading order in the long-wavelength perturbations, we have 

5.1 The Effective Fluid 

At super-Hubble scales, with k << H, the short-scale dynamics simply renormalize 
the background pressure, energy density, and anisotropic stress. We can determine 
the renormalization by evaluating [t^Ja as k ->■ 0. 

The background energy density and pressure are given by 

1 

Pes = lim - [r °] A Pes = lim - [r?] a (5.8) 

The background anisotropic stress is zero. 

5.2 Density and Velocity Perturbations 

To find the equations governing the dynamics of the effective fluid, we smooth the 
equations (I3.27H3.281 . Smoothing the continuity equation gives: 

pi + di(( Pl + Pt)vi) - vim + 3(H - $i)(pi + Pi) + 2v, • v,(p/ + P t ) 

= -diKps + P s )vi] A + [vidiP a ] A + 3[$ s (p s + P a )] A - [(p + P)d (v 2 )} s A (5.9) 



11 



From now on we will omit the "/" subscript from the variables to avoid cluttering the 
notation. Keeping only the terms linear in the perturbations 5i and vj, subtracting 
the homogeneous equation and dividing by p gives 

5= (l + w)(3$-V-v)-3H6(c 2 s -w) (5.10) 

To this order we can take w = jr. Similarly, smoothing Eq. fl3.28|) gives: 

v " * + (TTSj V ' v + V2 * = -W^) d ^ (511) 

where w = ■? as before and c 2 = j is the adiabatic sound speed. (The details of this 
calculation are somewhat involved and are therefore given in Appendix O) We can 
introduce the velocity potential v such that Vi = ikiV. 

cl , 1 / , „ o 9 



* + (if^ t '-* = -^)l c ^- 4 »«J (5 ' 12) 

Once again, our results differ from those given in [lj. We believe that the differences 
are due to extra simplifications that are made in (TJ. For example, Baumann et al. 
assume that the Newtonian approximation holds at scales k^ 1 << A -1 . 

The Einstein equations for the long-wavelength metric perturbations are given in 
Appendix IB. 21 



5.3 Scalar Field Perturbations 

Finally, we consider the long-wavelength scalar field perturbations by smoothing the 
Klein-Gordon equations. In order to simplify our calculations we will ignore the 
metric perturbations, which is a commonly used approximation when studying field 
dynamics during preheating [U [5j [Til 02] • Separating the fields 0, x m to homogeneous 
parts 0, x an d perturbations <50, 8x, and assuming that x — 0, we obtain the following 
equations for the perturbations: 

5(f) - V 2 (50) + 2H5<f> + a 2 m 2 5(f) + a 2 g 2 }5 X 2 (0 + 50) = (5.13) 
S'x - V 2 (5x) + 2H5 X + a 2 g 2 X ^ + 50) 2 = (5.14) 

Smoothing these equations and keeping only the terms linear in the long- wavelength 
perturbations gives: 

6'^i - V 2 (50O + 2K6'<j )l + aVSfa = -a 2 g 2 (<fe 2 ] A + [<^0 S ] A ) (5.15) 
S'xi - V 2 (5x0 + 2W X l + a 2 g 2 ^ 2 5xi = ~a 2 g 2 (20[5 Xs 50,] A + [5 X ^0 2 ] A ) (5.16) 
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We have obtained evolution equations for S(j)i and 5xi in terms of the 2 and 3-point 
correlation functions of the short-wavelength field perturbations. It would be inter- 
esting to find even approximate analytical forms for the correlation functions, so that 
the equations could be solved explicitly. 

6 Conclusion 

We have studied field dynamics during preheating in an unusual way, by focusing our 
attention on long-wavelength perturbations and integrating out the short-wavelength 
modes. As might be expected, we see that the homogeneous part of the inflaton 
field plays an important part in the evolution of both the <fi and x fields, as do the 
correlation functions between the perturbations 8(f) and 8\- We have also adapted 
and extended the formalism in [1] so that it can be applied to a universe whose 
matter content has pressure, and where the short-scale gravitational dynamics are 
not Newtonian. In contrast to the work in pQ, the only restriction we placed on our 
model was that the matter content should be modeled as a perfect fluid. 

There are several obvious ways in which this work could be extended. Firstly, 
of course, it would be desirable to find analytic solutions of the equations of motion 
given in this paper, and to compare them with the results of lattice simulations of field 
dynamics during preheating. We have also given the linearized Einstein equations, 
which govern the evolution of large-scale scalar metric perturbations during preheat- 
ing. Given that data collection is currently underway to find evidence distinguishing 
between various preheating models, it would be interesting to find ways of calculat- 
ing observable non-Gaussianities in the power spectrum of the metric perturbations 
using this formalism. Also, we have only considered scalar metric perturbations in 
this work: the next step might be to consider the evolution of large-scale vector and 
tensor perturbations. 

A Smoothing Bilinear and Trilinear Quantities 

Here we explicitly outline the calculations that lead to Eq. (14.511461) . Smoothing a 
bilinear quantity fg gives the same result as in [lj, but we find a different expression 
for [fgh]^. Therefore we will follow the approach given in [1] to derive [fg}\, before 
outlining our calculation for [fgh]\ and emphasizing the differences from the result 
in P. 
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Recall that we assumed the following useful properties for the Gaussian window 
function Wa(|x — x'|): 

^ A (x'-x) ? (x'-xy = ±5 ij (A.l) 

d f W A = -djW A = A 2 (x - x') J 'W A (A.2) 
d if d f WA = didjWA = -A 2 5ijW A + A 4 (x - x'^x - x')W A (A. 3) 

To smooth fg, we take its convolution with the window function: 

[fg] A = ^W A f(x')g(x'). (A.4) 

Splitting the fields /, g into long-wavelength modes fi, g\ and short-wavelength modes 
fa, 9s, we get: 

[fgh = [fwiU + [f s g s ]A + [fig.U + [f>9ih- ( A - 5 ) 

Since fi and gi are assumed to be small perturbations, and are long-scale, we can 
expand them in a Taylor series about x: 



/,(x') = /,(x) + <9^(x)(x' - x)< + -9,^/ z (x)(x' - x)*(x' - x)' + . . . (A.6) 
and similarly for gfj. This gives us 



[figiU = fm + ^ (yfi ■ + i//V 2 ^ + ^v 2 /*) + . . . (A.7) 

where the dots indicate higher derivative terms, suppressed by powers of k 2 /A 2 . To 
simplify the term [fig s } A , we first rewrite it as 

\fi9sU = Ma ~ [fWiU, (A.8) 

where the second term is given by (|A.7[) . To simplify the first term, we again use the 
Taylor expansion of fi, giving us: 

Ma = fm - difi[(x ~ x')Mx')]a + \mfi ■ t( x " x ') 4 (x - x')^(x')]a + . . . 

= fm + i hh ■ Vgi + lgiV 2 fl) +■■■ (A.9) 



A2 ^ Jl ^ 2 

Interchanging / and g in the above expression gives us \f s g\\ A . Thus we find that 

IfgU = fm + \f.g.]A + ^fi ■ v gi + . . . (A. io) 
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Smoothing a trilinear term fgh proceeds in the same way. We begin by splitting 
the fields f,g,h into long-wavelength modes fi,gi,hi and short- wavelength modes 
f„g„h a , giving: 

[fgh]A = [figihi] A + [f s g s h s ] A + [figh] A + [fgih) A + [fgh t ] A 

-IfgihiU-ifigh^-ifmhU (A.ll) 

Expanding fi,gi in Taylor series as before, we can smooth figih: 

[fWih] A = fmh - gAftK* - x') i /i(x / )] A - ffigMx - x / )^(x / )] A 
^^[(x-x^x-x'^a 

^gididjUix - x'^x - x')^] A + \did jgi [(x - x') ? (x - x') j '/i]a + • • • 

= figihi + gi — j- 2 — + fi — j- 2 — + hi — — 2 — 
V 2 f, V 2 q, 

+ gih^ + fih^ + . . . (A.12) 

We can also smooth figh: 

[figh] A = fi[gh] A - dMx - x') j gh] A + ~ " ~ x ')V]a + • • • 

f u^f\ h i f V9i-Vh t V/, -Vht Vfi ■ Vgi 
= figihi + fi[g s h s \ A + f — — 2 — + gi — — + h — — — 

+ 9ihl^ + £2 + —[g s h s ] A + ... 1 (A.13) 

where we have used flA.lOp . Substituting (IA.13]) and flA.12jl into f lA.lip . we find 

LfS^A = /tfi^i + L/sflAU + /j[sA]a + gi[f s h s ] A + hi[f s g s ) A 

+ jtifiVgi ■ vti, + • v/i, + fc,v/, ■ v#) 

+ ^i(V// ■ V[#A]a + Vgi ■ V[/A]a + Vh ■ V[f s g s \ A ) 

+ ^2^ 2 fi[9sh s ]A + V 2 </,[/A]a + V 2 hi[f s g s ) A ) + ... (A.14) 

The expression easily generalizes to smoothed polynomial quantities /1/2 • • • /n of any 
order n. 
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Baumann et al. claim in [TJ that smoothing the trilinear quantity pv l v J gives: 

[pvV} A = Pl v\v\ + [pv i y s ] A + pi — ^— l - + ... 

= PjU/Uj + [p^XU + [pj^XU + Pi — 1 + • • • 

= flt^ + [p^XU + PiK^k + ^ + Pi A 2 + • ' • 



(A.15) 



Comparing this to ( 1A.14j) . we see that our general expression for a smoothed trilinear 



quantity [fgh] A does not agree with Baumann et al.'s expression for [pv l v j ]a upon the 
substitution of p, v\ and t> J for /, and h, due to the presence of extra terms such as 
vl[p s vi}\ in our formula. As [fgh] A should be symmetrical in f,g, h, we believe that 
these terms should be present unless some additional assumptions or conditions are 
imposed. In our work, we will make no such additional assumptions. 

B The Einstein Equations 

B.l Velocity Expansion of the Einstein Equations 

As explained in Section 13.31 the energy-momentum pseudotensor t^ v is given by 

qNL 

> sr '-iS (B1) 

We expand these equations to order v 2 . The expansion of T^ v is given by Eq. fl3.191 - 
13.211) . As explained in Section [2, the metric perturbation $ is of order v 2 . However, 
each gradient V$ lowers the order in v by one, so terms of the form $V 2 $ are of 
order v 2 . We therefore work to first order in $, with the exception of such gradient 
terms. We end up with the equations 

V 2 $ - 3ft($ + WH) = -4vrGa 2 (r ° - T °) (B.2) 
di{® + m>) = knGa 2 Tl (B.3) 

<£ + H(2$ + if) + (H 2 + 2U)^ - ^V 2 ($ - *) = ^f^(r* - fi) (B.4) 

%TrGa 2 d l d j {T) - \5)tD (B.5) 



didj 



did^ - tt) - i%V 2 ($ - tt) 
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where the non-linear parts of the Einstein tensor are 

-a 2 (G° ) NL V$ ■ V$ + 4$V 2 $ (B.6) 

—a?(G i ) NL ~ (B.7) 
-a 2 {G)) NL ~ <Sj. V$ ■ V$ - 2V l $V i $ (B.8) 

to order v 2 . The homogeneous Einstein equations are 

u 2 = _^!f° = -*a°*- P (b.9) 

3 3 

U 2 + 2U = -^Tl = -^P (B.10) 
3 3 

B.2 The Einstein Equations in the Long- Wavelength Effec- 
tive Theory 

The long- wavelength metric perturbations are given by the Einstein equations, sourced 
by the effective energy-momentum pseudotensor. To linear order in the long- wavelength 
perturbations, we have 

V 2 $, - 3H($i + W&i) = -47rGa 2 (r ° - f Q °) = AnGa 2 5ip (B.ll) 
di($i + m>i) = 4ttGcl 2 t 1 = AnGa 2 (p + P)v\ (B.12) 

+ ft(2$, + %) + (H 2 + 2«)¥, - ^V 2 ($, - ¥,) = ^p-(rj - 7?) 



47rGa 2 2 

— - — c^pd/ (B.13) 



didj 



87rGa 2 ^(rj - i^r fe fe ) (B.14) 
3 



87rGa 2 (p + P)c4^ (B.15) 



where the components of r/f are 



r ° = - K]l - 87rGa2 (B.16) 

r* = fa + P,K + [fo + PM A (B.17) 

rj = ^ + 1(P + Pyvj]% - 1 k > kl 2 G J >- JA (B.18) 



Equations ( 1B.16tlB~T8|) differ slightly from those given in [TJ, as we have different ex- 



pressions for smoothed trilinear quantities. However, as the smoothed short- wavelength 



17 



perturbations can be written as a derivative expansion of long-wavelength variables, 
we find that to leading order in the perturbation variables, the Einstein equations are 
the same. 



C Smoothing the Momentum Conservation Equa- 
tions 

Here we show the explicit steps taken in smoothing the momentum conservation 
equation 

{g° v + u c u v )V^ = (p + PKVX + d a P + u a u^P = (C.l) 

which, when expanded to order v 2 , becomes: 

(p + P) (vi + v j d jVi + diV) + diP + VidoP + VitPdjP = (C.2) 

We will omit all the higher derivative corrections suppressed by powers of k 2 /A 2 that 
arise. Smoothing the term <9jP is easy: 

[d t p] A = m (c.3) 

Next we smooth the term (p + P)d^: 

[(p + P)d^] A = ( Pl + P l )d^ l + [(p s + P s )d^ s ] A (C.4) 
= ( Pl + P t )di% + ^-I-^V^Wa + [Psd&sU, (C.5) 

using the Poisson equation, V 2 ^ = A^Ga 2 p s . We then integrate by parts and use 
dj>W\ = —djW\ to get 

\{p + P)di*] A = (pi + Pi)d i V l + [Psd&sU + to^2 d A*,k*,k6i - 2*,&j] A (C.6) 
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Finally we smooth the remaining terms, integrating by parts and using the properties 
(1A.1HA.3) of the window function when necessary: 

[(p + P) (v { + V j d jVi ) + Vi P + VirfdjPU 

= d [{p + P)v% - J W A v\p + P) + J W A ((p + PydjVi + v'P + jvjdjP) 
= do[(p + P)v\ - j W A v l (-fy((p + Py) - 3(H - $)(p + P) - (p + P)(V)) 
+ J W A (p + P)v j d jVi 

= d [( P + p)v% + J w A (v i d 3 ,{{p + py) + ( P + pyd jVi ) 

+ I W A v* (3(H - $)(p + P) + (p + P)(v 2 ) 



(C.7) 



d [{ P + py] A + d 3 J w a ( p + p)vV 

+ / W A v* (?>{% - $)(p + P) + (p + P)(t;' 2 ; 



Using the smoothed continuity equation ( 15. 9ft . we can substitute for p + P, which 
gives: 

[(p + p) + + ^p + w t va,-p] A 
= ( Pl + p t y + d [( Ps + p s y s ] A + o 3 J w A ( P + p)vv 

+ «/ (-dAiPi + PiU) - 9A(Ps + PsKU + Pi + v\d 3 P { + [vid 3 P s } A ) (C.8) 
Adding Eg. (10311081) gives: 

( Pl + p t y + ( Pl + Piydjvi + v\P x + vjvjdjPi + ( Pl + p l )d i ^> l + m 



>» ^ j 



87rGa 2 



+ ^-[(p s + P s K] A -^[^P 



a A 



J A 



- OoKps + PsH) A - d 3 [{ P + pyyy A - [pa* 



s A 



(C.9) 



Keeping only the leading terms linear in the matter and velocity perturbations, we 
get: 



P BP 
(1 + wy + v\-± + (1 + w)di% + ^ = d 3 [r l3 } s A 
Pi P 



(CIO) 
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